Distinguishing Modified Gravity from Dark Energy 
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The acceleration of the universe can be explained either through dark energy or through the 
modification of gravity on large scales. In this paper we investigate modified gravity models and 
compare their observable predictions with dark energy models. Modifications of general relativity 
are expected to be scale-independent on super-horizon scales and scale-dependent on sub-horizon 
scales. For scale-independent modifications, utilizing the conservation of the curvature scalar and a 
parameterized post-Newtonian formulation of cosmological perturbations, we derive results for large 
scale structure growth, weak gravitational lensing, and cosmic microwave background anisotropy. 
For scale-dependent modifications, inspired by recent f{R) theories we introduce a parameterization 
for the gravitational coupling G and the post-Newtonian parameter 7. These parameterizations 
provide a convenient formalism for testing general relativity. However, we find that if dark energy 
is generalized to include both entropy and shear stress perturbations, and the dynamics of dark 
energy is unknown a priori, then modified gravity cannot in general be distinguished from dark 
energy using cosmological linear perturbations. 



I. INTRODUCTION 

Cosmic acceleration has been revealed by measure- 
ments of the redshift-distance relation x(^) where x is the 
comoving radial distance. The Hubble expansion rate fol- 
lows from H{z) = {dx/dz)~^ (in units where c = 1). This 
determination assumes only that the observable universe 
is adequately described by the Robertson- Walker metric, 
an assertion that is testable empirically ^, ^ and does 
not imply the validity of general relativity (GR). 

The inference of dark energy follows once the Ein- 
stein field equations of general relativity are imposed 
on the Robertson- Walker metric yielding the Friedmann 
equations. These equations imply that a stress-energy- 
momentum component with negative pressure is needed 
to explain cosmic acceleration. This substance may be 
vacuum energy (i.e., a cosmological constant, giving rise 
to the ACDM model) or a scalar field 0]. The dark en- 
ergy equation of state for uniform expansion, p{p), can 
be determined from measurement of p{z), which itself fol- 
lows from H(z) combined with the first Friedmann equa- 
tion. Measuring w{z)= p/p is the primary goal of dark 
energy experiments |4'|. 

Another possibility is that general relativity requires 
modification on large distance scales and at late times in 
the universe. In this case cosmic acceleration would arise 
not from dark energy as a substance but rather from 
the dynamics of modified gravity. Modified gravity is 
not particularly attractive theoretically, but the observed 
cosmic acceleration is so surprising that all plausible ex- 
planations should be considered. 

Observations of the cosmic expansion history cannot 
distinguish dark energy from modified gravity Q. 
Testing gravity requires exploring the evolution of spa- 
tial inhomogeneity (e.g. [J Q and references therein). 
Modified gravity theories must pass tests within the so- 
lar system and in relativistic binaries Q. They are ex- 
pected to show significant departures from general rela- 
tivity only on cosmological distance scales. The combi- 
nation of cosmic microwave background anisotropy, weak 



gravitational lensing, and the growth of clustering of dark 
matter and galaxies provides an opportunity to discrim- 
inate between dark energy and modified gravity. 

Performing cosmological tests of modified gravity re- 
quires a set of predictions. There are two approaches to 
generating these predictions. In the first, a theory is spec- 
ified by its Lagrangian (or other fundamental description) 
which provides the equations of motion for both homoge- 
neous expansion and cosmological perturbations. A class 
of theories can be specified by giving a Lagrangian with 
free parameters, e.g. f{R) theories where R is the Ricci 
scalar @. 

A second approach is inspired by the parameterized 
post-Newtonian framework for solar system tests jTo| . 
Here one begins with the solution of the gravitational 
field equations (i.e., the metric) instead of the La- 
grangian. Several authors have recently adopted this 
framework or a similar one [ll|, [l^, [H, [3, [3] ■ The dif- 
ficulty here is to find a good "Newtonian" description in 
cosmology to which one adds "Post-Newtonian" parame- 
ters. Metric perturbations in the scalar sector governing 
the growth of cosmic structure are characterized by two 
spatial scalar fields, $ (the Newtonian potential) and "if 
(the spatial curvature potential). Even if we introduce 
the relationship ^! — 7$ with Eddington parameter 7 
[l^ . there remains one unknown function of space and 
time. In the solar system case, by contrast, $ = —GM/r 
is known to provide an excellent approximation to plan- 
etary dynamics. 

On solar system scales, and even within galaxies, one 
can use test-particle orbits to determine $, and light rays 
to determine $ -I- 5*. This comparison yields impressive 
limits on I7 — 1| in the solar system Q. However, on a 
scale of several kpc, gravitational lensing combined with 
stellar dynamics in elliptical galaxies yields a current best 
result I7-II = 0.02±0.07 [l^. At the scale of Gpc where 
dark energy appears to drive accelerated expansion, there 
are no longer any bound test-particle orbits to measure 
$, so a different approach, based on cosmological pertur- 
bation theory, is needed. 
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Previous work in the cosmological parameterized post- 
Newtonian framework has either assumed that some of 
the Einstein field equations remain valid with modified 
gravity [l^ or has examined the dynamics of individual 
theories, e.g. [l^. Neither approach is ideal. One would 
prefer to sample all possible theories in a broad class, and 
for each theory to constrain the potentials by a consis- 
tency condition that does not assume general relativity 
or any particular modification thereof. 

Such a consistency condition was found recently in 
Ref. [19| for the long-wavelength perturbations of a 
Robertson- Walker spacetime. This result was derived 
assuming that gravity is described by a classical four- 
dimcnsional metric theory having a well-defined infrared 
limit (i.e., the theory is well behaved for very long wave- 
length perturbations). For practical application, assump- 
tions must also be made about the background spatial 
curvature and entropy perturbations. Assuming an in- 
flationary or equivalent origin of perturbations, long- 
wavelength isentropic perturbations are imprinted in the 
spatial curvature on a flat background. In general rel- 
ativity, these spatial curvature fluctuations, represented 
by the gauge-invariant C field of Bardeen et al. [2(| or 
the TZ field of Lyth 2l|, are time-independent in the 
long- wavelength limit. Ref. [l^ presented a derivation of 
the conserved curvature perturbation (calling it k) mak- 
ing no assumption about the field equations except that 
they have a well-defined infrared limit. Physically this 
means that curvature perturbations are small and that 
all waves propagate causally. In what follows, the curva- 
ture perturbation introduced in Ref. [l^ will be called 
C although its definition differs from that of Ref. [20l |. 
For long wavelength perturbations on a flat background, 

n = c. 

In the long-wavelength limit all cosmological pertur- 
bations factorize into functions of time multiplying the 
curvature perturbation or its spatial derivatives. This is 
true in GR and in modified gravity theories that are well- 
behaved in the infrared limit. One might naively expect 
this factorization to hold only on scales larger than the 
Hubble length. In general relativity, however, signals in 
the scalar sector propagate at the speed of sound, not the 
speed of light 2z\ , leading to conservation of C on scales 
larger than the Jeans length. 

To satisfy solar system tests, modified gravity theories 
for cosmic acceleration must introduce a length scale Lq 
below which general relativity is recovered. This length 
scale might be associated, for example, with the dynam- 
ics of new scalar degrees of freedom. The value of this 
length scale, compared with the size of the systems inves- 
tigated, plays a crucial role in characterizing the behavior 
of modified gravity theories. 

If Lg is much smaller than the length scales over which 
linear cosmological structure formation is measured (e.g., 
Lq = 1 Mpc), then the factorization of cosmological per- 
turbations on scales larger than the Jeans length remains 
valid. We denote this case scale-independent modified 
gravity. These theories are like GR in that the curvature 



perturbation is conserved for the relevant length scales. 
This condition yields a great simplification of the dynam- 
ics, reducing cosmological perturbations to quadratures. 

In GR, waves propagating in the scalar sector travel 
only at the speed of sound, so that scale-dependence 
of transfer functions arises only below the Jeans length. 
Modified gravity theories, however, typically have addi- 
tional fields supporting waves that travel at the speed of 
light. In this case, Lq is the Hubble length and the fac- 
torization of cosmological perturbations no longer holds. 
Theories of this type are called scale-dependent modified 
gravity theories. Now two quantities, 7 and the gravita- 
tional coupling (the generalization of Newton's con- 
stant in the Poisson equation for are needed to char- 
acterize gravity, and both will vary with length scale as 
well as with time. Even such complicated models can still 
be approximated by parameterizations, as we will discuss 
below. Ref. [TTl found a way of bridging super- and sub- 
horizon modifications to GR. Our parametrization, while 
not as general, is simpler because it only involves a few 
free parameters and no free functions. 

Because we do not start with a Lagrangian, we cannot 
explain cosmic acceleration. We take the cosmic expan- 
sion history as given from observations. Rather than pro- 
viding a complete theory of modified gravity, we provide 
a framework for observational tests of gravity in cosmol- 
ogy. 

This paper is organized as follows. Section [III describes 
the curvature perturbation and its use to build scale- 
independent modified gravity theories. Section Hill works 
out the growth of structure on sub-horizon scales, shows 
that the Poisson equation is modified, and derives results 
for cosmic microwave background anisotropy and weak 
gravitational lensing for scale- independent modified grav- 
ity theories. Section IIVI then examines the sub-horizon 
behavior of a currently popular class of theories known 
as f{R) models and shows that they are scale-dependent. 
Section lVTl considers the alternative hypothesis that dark 
energy has a peculiar stress tensor while gravity is gov- 
erned by GR. Finally, results are summarized and con- 
clusions are presented in Section IVIII 



II. GRAVITY AT LONG WAVELENGTHS 

Our starting point is the perturbed Robertson- Walker 
metric in conformal Newtonian gauge p3| : 



a^{t)[-{l + 2^)di^ + (1 - 2^)"f^jdx' dx' 



(1) 



where t is conformal time, a{t) = 1/(1 + z) is the expan- 
sion scale factor, and "fij (x, K) is the three-metric for a 
space of constant spatial curvature K, e.g. ^ijdx'^dx^ — 
dx^ _j_ r^(x,i^)dri^ where r(x^K)^/ll = sin(x-\/^) for 
K > {) and is analytically continued for K < Q. Note 
that different conventions appear in the literature for the 
metric perturbations: $ = V^hu = and ^ = — ^Hu = 4> 
where (^Hu,'i'Hu) are the potentials of Ref. [llj] and 
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(■0, <f) are the potentials of Ref. [23|. Linear perturbation 
theory is assumed to be valid throughout this paper. 

The evolution of the scale factor can depend, in prin- 
ciple, on any quantities characterizing the geometry and 
composition of the Robertson- Walker background, for ex- 
ample the spatial curvature K and the entropy density 
(or equivalently, parameters characterizing the equation 
of state). We neglect entropy perturbations and con- 
sider only curvature perturbations on a flat {K ~ 0) 
background. Assuming that the unknown gravity the- 
ory has a well-defined infrared limit obeying causality, 
long wavelength curvature perturbations should evolve 
like separate Robertson- Walker universes. In this case 
it is possible to transform to a new set of coordinates, 
t ^ t — a{t) and x ^ x(l + C) where C is constant 
and d = $ -I- — such that the new line element is 
eq. H]) with $ = = and having spatial curvature 
K{1 + 2C). Thus, C is one-half the spatial curvature per- 
turbation. Enforcing the coordinate transformation leads 
to the consistency condition [l^ 

(2) 

where TL = a/ a = aH and k is the comoving wavenum- 
ber. Although Ref. [Si states that large-scale shear 
stress is neglected in this result, in fact eq. ([2]) is valid 
for fc — > in general relativity (and presumably in modi- 
fied gravity theories) even if shear stress is present. The 
curvature term KjH? has been computed assuming the 
Friedmann equation is valid; in modified gravity theo- 
ries this term might be different but it must vanish when 
K = Hereafter we assume if = and drop the curva- 
ture term. 

Eq. ([2|) may be regarded as a definition of the curva- 
ture perturbation C for arbitrary theories of gravity. For 
long wavelengths, C is independent of time. Sound waves 
in the matter sector or wave propagation in the modi- 
fied gravity sector cause C to change with time on small 
scales. These changes are implied by the neglected terms 
proportional to fc^C iii eq. ([2]). For now we ignore such 
terms, in effect assuming that both the Jeans length and 
Lc are smaller than the cosmological scales of interest. 

During the radiation-dominated era the Jeans length 
is comparable to the Hubble length, and C, (and <i> and ^') 
is damped for scales smaller than the Jeans length. We 
assume that during this early period of evolution general 
relativity is an excellent approximation so that the damp- 
ing is well described by the transfer functions computed 
using standard codes 0, HI]. When modified gravity 
becomes important at low redshift, the Jeans length has 
dropped to a few Mpc or less. In practice, we modify 
CMBFAST only for z < 30 and then do so in such a way 
as to enforce eq. ^ with corrected from its primeval 
value using the GR transfer function at z = 30. 

Assuming a well-defined infrared limit, the time and 
space dependence of perturbations must factorize for 



wavelengths longer than the Jeans length or Lq^ e.g. 

<I>(k,t)=i^(a)C(k)+0(fc2C) (3) 

where k is the wavevector. Factorization implies that the 
ratio of the two gravitational potentials depends only on 
time as fc ^ 0. Therefore we may write, for any causal 
theory of gravity having a well-defined infrared limit, 

vl/(k,i) =7(a)$(k,t)-f 0(fc20 . (4) 

In modified gravity theories, 7(a) is the only degree of 
freedom important for long-wavelength scalar perturba- 
tions. The conditions of causality and a well-defined 
infrared limit greatly restrict the dynamics of modified 
gravity theories. 

We now make the key assumption that the terms pro- 
portional to fc^C in eqs. ©-(H]) can be neglected not only 
on super-horizon scales k < H^^ but also, as they can be 
in general relativity, on sub-horizon scales down to the 
Jeans length. This assumption defines a class of theories 
we call scale- independent modified gravity models. 

Under these assumptions, modified gravity is com- 
pletely specified on large scales by the scale-independent 
function 7(a) . At high redshift when dark energy is unim- 
portant we require 7 — > 1 in order to retain the success 
of general relativity in explaining the cosmic microwave 
background anisotropy [26| . Thus we adopt the following 
parameterization for scale-independent modified gravity: 

7(a) ^1 + I3a% (5) 

where (3 and s > are constants. Eqs. (H))^® now give 

7F(a) = 0-2^7^'/^) r 07-^'/^) ^(-]da. (6) 
Jq da \T-LJ 

Changing the lower limit of integration introduces a 
rapidly decaying solution which we ignore. 

In general relativity with negligible shear stress, 7 = 
1. When a component with constant equation of state 
parameter w > — i is dominant, a oc with n = 2/(1 -|- 
iw) yielding F — {i + 3w)/(5 + 3w). Thus, for long 
wavelengths the potential drops from $ = |^ during the 
radiation-dominated era to = |^ during the matter- 
dominated era. 

Of greater interest here is the evolution of the po- 
tentials during the matter-dominated era with modified 
gravity parameterized by ([5]). Figure [1] shows the results 
for s = 1 and s = 3 as well as the GR case /3 = 0. 
The background expansion history is chosen to match 
GR with Vljn = 0.284 and a cosmological constant with 
51a = 0.716. The choice s = 3 matches Caldwell et al. in 
the limit j3 [12(]. However, as we will see below, our 
results differ from theirs because they assumed the valid- 
ity of some components of the Einstein field equations, 
while we instead required consistent causal evolution on 
large scales. This difference will be discussed further be- 
low in Section IVll 
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III. OBSERVABLES FOR 
SCALE-INDEPENDENT MODIFIED GRAVITY 

With the time evolution of the metric in hand for 
long wavelengths we are now able to calculate observable 
quantities for scale-independent modified gravity theories 
parameterized by the constants (/?, s). The effects con- 
sidered here are the growth of structure in the dark mat- 
ter, microwave background anisotropy, and weak gravi- 
tational lensing. 



A. Growth of structure 

Until now, no assumptions have been made about dy- 
namics in the matter sector except for causality and 
consistency with a spatially homogeneous, uniformly 
expanding Robertson- Walker solution. To follow the 
growth of structure we must specify how the matter 
fields are coupled to gravity. Here we assume that the 
dark matter obeys the weak equivalence principle, i.e. 
collisionless dark matter particles follow geodesies. This 
choice explicitly forces scalar-tensor theories to the Jor- 
dan frame in which matter fields are minimally coupled 
to gravity. 

In the conformal Newtonian gauge, on sub-horizon 
scales where |i5| ^ \^\ with i5 = bpjpra, and pm is the 
average mass density, cold dark matter fluctuations obey 
the evolution equation 



(5 + H<5 = -k^<^ 



(7) 



This equation follows from particle number conserva- 
tion and geodesic motion or, equivalently, from energy- 
momentum conservation. The density perturbation field 
can be written as 



FIG. 1: Evolution of the scalar potentials $ and ^ in the 
long wavelength, scale-independent limit, assuming a. C, = \ 
normalization. Modified gravity effects, parameterized by eq. 
((S)), arise later for larger s. The GR model assumes 7 = 1 and 
a cosmological constant. For 7 < 1 the Newtonian potential 
grows, unlike general relativity with a cosmological constant. 



5(k,i) = -A:2i?(a,fc)G(k) 



(8) 



where is the curvature perturbation at a = which 
we take to be 2; = 30 so that the Jeans length is smaller 
than the scales of interest and modified gravity has not 
yet become important. For a > a^, $ factorizes and 
eq. ([7]) can be reduced to quadratures for D{a,k). With 
initial conditions D{a, k) = Di{k) and dtD{a, k) = Di{k) 
at a = a, , the solution is 



Figure [T] shows that the Newtonian potential $ is en- 
hanced and the spatial curvature ^' = 7$ is diminished 
for 7 < 1, compared with general relativity. For s — 3 
the modifications occur later because I7 — 1| is smaller 
at earlier times. The quantitative results depend on the 
validity of eq. ^ but this qualitative behavior (the New- 
tonian potential being enhanced for 7 < 1) should persist 
in general. 



where 



D(a, k) = Dia) + D,{k) + a,y{a)D,{k) 



D{a) = y / j^da- I ^ da , 



via) 



° da 



(9) 

(10a) 
(10b) 



The function y{a) asymptotes to a constant but D{a) cx a 
for a ^ Oi in the matter-dominated era. Thus the late- 
time solution for density perturbation growth factorizes, 
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FIG. 2: Evolution of the logarithmic density perturbation 
growth rate d\nD/d\na. Models with 7 < 1 have enhanced 
growth relative to the ACDM (GR) model. 



D{a,k) — D{a) in cq. ([H]). This perturbation growth is 
often represented as a function of redshift by defining 
g{z) = D{a)/D{l) with a = 1/(1 + z). 

Figure shows the logarithmic derivative of the den- 
sity perturbation growth versus time for our parame- 
terized modified gravity models. In the ACDM model, 
dlnD/dlna w ['[}m{a)f^^^ \2^- The transition to a cos- 
mological constant-dominated universe leads to a sup- 
pression of growth. If instead gravity is modified, the 
growth rate can be increased or decreased relative to the 
GR case. The qualitative effects are easy to understand. 
We already saw that models with 7 < 1 (/3 < 0) get an 
enhanced Newtonian potential A stronger potential 
increases the gravitational force on density perturbations 
leading to an enhanced growth rate. 

The simplest test of growth of perturbations is the total 
perturbation growth by redshift zero, which is character- 
ized by the variance of density fluctuations in spheres of 
radius Rs ~ 8 Mpc, 



(11) 



where W{x) — 3ji{x)/x. The power spectrum of matter 
density fluctuations is 



D2(z = 30) ' 
(12) 

where A-ji is the amplitude of the initial scalar curvature 
fluctuations on scale fco and Tm is the transfer function for 
matter fluctuations in the synchronous gauge relative to 
( = 71 computed by CMBFAST (which accounts for the 
suppression of growth during the radiation-dominated 
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FIG. 3: Contour plot of crs (/3, s)/cr8 (ACDM) with /J running 
along the horizontal axis and s running along the vertical axis. 



era). We adopt A^, = 2.4 x 10"^, fco = 0.002 Mpc"\ 
and Us = 0.958 26]. 

Our modification of gravity is scale-invariant in that 
the fc-dependence of the dark matter power spectrum is 
unchanged relative to GR. Thus the amplitude of den- 
sity perturbations as depends on modified gravity only 
through the enhancement or diminution of growth shown 
in Fig. H 

The specific results obtained here assume that the 
gravitational potentials factorize on scales larger than a 
few Mpc. If this is true, then CMB and galaxy clustering 
results on all scales can be fit by a single modified gravity 
theory with one function 7(a) or, equivalently, D{a). If, 
however, modified gravity introduces a length scale Lq 
intermediate between i?g and the Hubble length, then 
the growth of structure will depend on wavenumber in a 
way not described by a scale-independent D{a). In Sec- 
tion |V] below we consider an alternative scale-dependent 
parameterization of modified gravity and examine the ob- 
servable consequences. For now, we assume Lq < Rs or 
equivalently that super-horizon relations apply, as they 
do in GR, to sub-horizon scales all the way down to the 
larger of the Jeans and nonlinear lengths. 

Figure [3] shows a contour plot of erg (normalized to the 
GR case) for different choices of the modified gravity pa- 
rameters. As expected, smaller values of 7 (i.e., /3 < 0) 
lead to larger amplitude. Thus, modified gravity changes 
the amplitude of galaxy clustering relative to the CMB, 
and could explain any apparent discrepancy between the 
values of erg inferred from CMB analysis and galaxy clus- 
tering or lensing measurements. 
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B. Modified Poisson Equation 

Rearranging eqs. Q and ([8]), we arrive at a modified 
Poisson equation relating the Newtonian potential to the 
overdensity: 



D(a) 



(13) 



The space curvature potential obeys a similar Poisson 
equation, 



(14) 



where G,f (a) = 76$ (a). Plots of G$ and are shown 
in Figure 3) Their time dependence is dominated by the 
potentials F{a) and 7(a)F(a) since 6 is less sensitive to 
our modified gravity parameters (/?, s). As a result, we 
see the same qualitative behavior as in Fig. [T] Models 
with 7 < 1 produces a greater value of G$ relative to the 
ACDM model, and larger values of s produce larger late 
time behavior. 

In GR, the gravitational coupling G is constant. In 
many alternative theories of gravity, the strength of grav- 
ity varies with time (and also with place, for length scales 
smaller than Lq). Time- varying G is well known for 
scalar-tensor theories [2^, [2^, but we find that it is a 
generic feature of all modified gravity theories with 7 7^ 1 
on cosmological scales. 

The time- variation of G, represented by G/G, has been 
severely constrained by measurements in the solar sys- 
tem, in stars, and in the early universe [sot - Limits on 
larg er scales are provided by the microwave background 
|3lj . The CMB acoustic peak structure will be unaf- 
fected if variations occur only long after recombination. 
Modified gravity explanations of cosmic acceleration sug- 
gest the need to constrain G/G on large length scales in 
the late universe. It would be very interesting to know, 
for example, to what extent the structure of clusters of 
galaxies and their X-ray emission can be used to con- 
strain G/G. 

Note that the method used to derive our modified Pois- 
son equations is roundabout. Had we started with a 
Lagrangian, the gravitational field equations would di- 
rectly yield the gravitational coupling strength. Because 
we started with a phenomenological description of modi- 
fied gravity, we instead deduce the dynamics of this cou- 
pling from the requirements of causal evolution and the 
weak equivalence principle. On small scales our treat- 
ment breaks down as modifications of general relativity 
must become scale-dependent. Nonetheless, the motiva- 
tion to investigate limits on G/G on scales much larger 
than the solar system remains valid. 



C. CMB temperature anisotropy 

Modified gravity (or dark energy) affect the microwave 
background only at late times through the integrated 
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by the potentials and exhibits the same features shown in Fig. 
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Sachs- Wolfe (ISW) effect: 
AT, 



(n) = /(«' + $) dx 



(15) 



where x is the comoving radius and h is the photon di- 
rection. The ISW effect arises when the gravitational po- 
tentials change with time, as occurs during transitional 
periods in cosmic evolution. One such contribution oc- 
curs during the transition from radiation to matter domi- 
nation. The other contribution is occurring today during 
the current transition to an accelerating expansion. The 
physics governing the matter-radiation transition is well 
explained by GR, while the physics governing the transi- 
tion today is (for the models considered here) dependent 
on modified gravity parameters. Because the recent ISW 
effect arises relatively nearby (z < 2), it shows up only 
at large angular scales. The ISW contribution from some 
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FIG. 5: The temperature anisotropy power spectrum for dif- 
ferent parameter choices of modified gravity. 



FIG. 6; Ratios of weak lensing shear correlation C+, to its 
value for GR, for several sets of modified gravity parameters. 



modified gravity models has been studied in several re- 
cent papers [gI. H^. 

We computed the CMB temperature anisotropy spec- 
trum by modifying CMBFAST [H to replace the ISW 
contribution of the ACDM model with that for our mod- 
els assuming the factorization of the potential. The 
results are shown in Figure [5l As expected, only the 
low-order multipoles are affected. Models with higher s 
produce larger changes because they lead to larger time 
derivatives. Models with 0.2 < 7 < 1 (-0.8 < /? < 0) 
produce less anisotropy because of destructive interfer- 
ence between the ISW and primary anisotropy contri- 
butions. Although decreasing the quadrupole moment 
improves agreement with observations, little statistical 
weight can be given to this conclusion because the mod- 
ifications, at least for 0.2 < 7 < 1.5, are smaller than 
cosmic variance. However, cross correlating the CMB 
with galaxy surveys could potentially be a more discrim- 
inating probe of modified cavity [12, 

It was recently found [34l| that our results are consis- 
tent with recipe Rl of Caldwell et al. However, we 
expect differences from our Figure [S] for recipe R3 of [12 
since C is not conserved on large scales in this scheme. 
This difference will be discussed below in Section IVTl 



D. Weak lensing 

Metric perturbations <I> -I- affect both the energy of 
photons (ISW effect) and their direction of travel (gravi- 
tational lensing). Gravitational lensing causes both mag- 
nification (or de-magnification) and differential stretch- 
ing (shear) of background images. The correlation func- 
tion or power spectrum of weak gravitational lens shear 
is an observable measure of large-scale structure. The 



weak lensing power spectrum is given by (3£ 



TDK 



X 



I 

X 



X) , (16) 



where 



Pv 



(1+7)^ 



and 



fc3 



W{x) = 



F^z) 



dx' 



X - X 
X' 



30) ' 
(17) 



(18) 



Here Xoo is the comoving distance to 2: = 10 (the results 
change by a negligible amount if the maximum redshift 
lies anywhere between 6 < 2: < 15), T$ is the transfer 
function of the Newtonian potential relative to C com- 
puted at 2; = 30 using CMBFAST and r]{x) is the radial 
distribution of sources, normalized with J 'r]{x) dx — 1- 
Note that these formulae assume a flat space. Our lensing 
analysis uses the source distribution 



T]{z) cx z^exp[-(1.41z/zmcd) 



1.5i 



(19) 



with Zmed = 1-26 [361 . This distribution approximates 
the galaxy redshift distribution of the COSMOS survey, 
if there were no clumping. 

Measurements of weak gravitational lens shear, for 
galaxies separated by angle 9 on the sky, provide an es- 
timate of shear correlation functions including 



dl . 



(20) 



Modifying gravity changes F{z) thereby changing C+{9). 
Figure [H] plots C-f (0) for modified gravity, normalized at 
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each 9 by its GR value. At z = 1, 1 Mpc corresponds 
to 2.15 arcmin. Hence, for some of the scales shown m 
Figure [6] structures are nonlinear and thus beyond the 
regime of validity of the current framework. However, 
a scheme that takes a linear power spectrum to a non- 
linear power spectrum would correct this flaw [s^l- On 
such small scales, our assumption of scale-independent 
modified gravity may also be invalid, so Figure [H] should 
be regarded as suggestive, but not definitive, of modified 
gravity effects on weak lensing. 

As expected, models with 7 < 1 have larger shear cor- 
relations because they have a larger F{z) and therefore 
more growth of structure (despite having a smaller I-I-7). 
For angular scales less than about 10 arcmin, the effect 
is almost equivalent to a constant change in the normal- 
ization of the power spectrum, i.e., in the value of ug- At 
larger angular scales, the redshift-dependence of F{z) at 
small redshift translates to a dependence on distance and 
hence on angular scale, however this is in a regime where 
the shear correlations are small and difficult to measure. 
Thus, scale-independent modified gravity theories pre- 
dict an amplitude of weak lensing different from GR with 
the same CMB primary anisotropy. The acoustic peak 
amplitudes tightly constrain In principle, measure- 
ments of (78 based on the CMB acoustic peaks (which are 
unaffected by modified gravity) could differ both from 
measurements based on galaxy clustering [which depend 
on D{z)] and those based on weak lensing [which depend 
on F{z)\. Current error bars are inconclusive [2^ but 
this comparison of different ag values could eventually 
provide a powerful test of GR. 



IV. COMPARISON WITH f{R) THEORIES 

Substantial work has already been done investigating 
modified gravity effects for f{R) theories. Here we con- 
sider such theories where the Ricci scalar R is replaced 
by i? + /(i?) in the Einstein-Hilbert action, and where 
the action is extremized with respect to the metric. In 
these models, the field equations are generically fourth- 
order. In effect, modified gravity introduces a new prop- 
agating scalar degree of freedom coupled to gravity, the 
scalaron Jr = df /dR [llj. The Compton wavelength of 
the scalaron imprints a physical length scale, which is 
made dimensionless by combining with the wavenumber 
k: 



Q 



3P I 



B.R 



(21) 



where fan = d^f/dR^. Several papers have recently 
discussed cosmological perturbation evolution for metric 
f{R) theories 39, 40, 41]; our notation most closely fol- 
lows that of ref. [41[ except our potentials $ and ^' are 
exchanged from theirs. 

In f{R) theories, C is conserved on super-horizon scales 
[ll|. However, the scalaron obeys a nonlinear Klein- 
Gordon equation with two length scales: the Hubble 



length and the scalaron Compton wavelength. The in- 
teresting case for large-scale structure is the quasi-static 
regime of linear, sub-horizon perturbations (fc^ 3> Ti.^) 
where 41 1 



l + V3 + 3Q 



5 , 7 



3 + 2Q 
3 + 4Q 



(22) 



Differentiating eq. ([2]) and substituting eqs. (O and ([22]) 
along with the background evolution equations (5) and 
(6) of ref. [41] for a universe containing only nonrela- 
tivistic matter and (optionally) a cosmological constant 
yields 



where 



U = 



aHd_ f V^B 



2QB 
3 + 2Q 



(23) 



(24) 



and 



V 



d_ 

di 



TB d 



4 © 



where we have defined the auxiliary variables 



F 

B 



1 



1 



\ d 




di 





-1 



3-h3Q 



2(1 + fR)n^ 



(25) 

(26a) 
.(26b) 



General relativity with a cosmological constant corre- 
sponds to the case / = 2A, Jr = 0, and 7 = F = 1, yield- 
ing U = V = Q. Thus, C, is conserved even on sub-horizon 
scales in a ACDM universe. However, this is no longer 
true if /fl ^ 0. Two distinct effects modify the curvature 
perturbation. First, the 1 + Jr factor in (|^^ modifies the 
evolution on sub-horizon scales. In practice, this effect is 
small if ^ 1, as is favored by galactic structure con- 
siderations [40]. In this case, the background expansion 
history is nearly identical to GR with a cosmological con- 
stant and gravity is significantly modified only at wave- 
lengths approaching the scalaron Compton wavelength, 
where Q ^ 1. For long wavelengths such that Q ^ 1, 
7 « 1 — |Q and the corrections introduced to eq. ([2|)-(j4|) 
by scalaron dynamics are O(fc^). The treatment given 
in the preceding sections remains valid for ]/fl] <C 1 and 
Q <C 1. However, this limit corresponds to general rela- 
tivity. Unfortunately, the treatment presented in Section 
ITTl which was based on a scale-invariant modification of 
gravity, breaks down just where f{R) theories begin to 
deviate significantly from GR. 

The f{R) models generically have 7 — 1 oc A:^ for sub- 
horizon wavelengths longer than the scalaron Compton 
wavelength. These models have a scale-dependent modi- 
fication of gravity. Although we cannot use the results of 
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Section |TT] to describe them, it is still possible to param- 
eterize scale-dependent modified gravity models so as to 
obtain useful results for the sub-horizon growth of large- 
scale structure. A simple parameterization inspired by 
f{R) theories is presented in the next section. 



V. SCALE-DEPENDENT MODIFIED GRAVITY 

For a wide class of theories, modified gravity leads 
generically to a Poisson equation with variable gravi- 
tational coupling. In the scale-invariant modifications 
of Section [ITl the Newton constant is replaced by the 
time- varying (t) which follows from the scale-invariant 
potential ratio 7(i). In scale-dependent modified grav- 
ity theories, on the other hand, G^{k,t) and "f{k,t) — 
Gip / G$ are functions of length scale as well as time, and 
there is no simple relation between them. Thus, more 
parameters are needed to characterize such theories [3] . 

Despite their generality, f{R) theories with /r 1 
have, for a wide range of sub-horizon length scales, a 
very simple form for and 7 given by eq. (P^ . To 
arrive at a simple phenomenological model we simplify 
the time dependence as follows: 



G$ l + Qfifc^a" 
G ~ 1 -I- a2k'^a'' 



1 



(27) 



We assume that these relationships hold only in the linear 
regime of cosmological density perturbations, and that 
G$/G — + 1 and 7 ^ 1 on solar system scales. We also 
require GR to hold at early times, implying s > 0. Eq. 
(I?7)) describes f{R) theories with ^ 1 if ai = = 
2/3i = — As a simple post-Newtonian 

model we will now assume that {ai, a2, P2) are arbi- 
trary constants with units of length squared. In order to 
ensure that G.^/G and 7 are finite for all fc, we require a2 
and P2 to be non-negative. Moreover, we need G* > 
in order to ensure that gravity is attractive. Hence, ai 
must be non-negative as well. This scale-dependent pa- 
rameterization has a different dependence on length scale 
than that of Amin et al. [l^. It is chosen to reproduce 
the scale-dependence of f{R) theories. For some mod- 
ified gravity theories 7 = 1, e.g., Einstein plus Yukawa 
gravity. For this model G $ / G in eq. ([77)1 is multiplied by 
an overall factor 0:2/0:1 [la | so that the deviation from 
Einstein gravity shows up only at large distances. 

The class of theories considered here has at least three 
physical length scales: the horizon scale a/Ti, the transi- 
tion scale a-'^+^/^Y^oT where gravity changes strength (for 
simplicity, we consider models where the and (3i are all 
of comparable magnitude) , and the nonlinear length scale 
for structure formation (e.g., approximately 10 Mpc to- 



day). If 



s/2 



'Oi and a^^^^'^y/Pi are smaller than the 



nonlinear scale, then for purposes of large-scale struc- 
ture formation, gravity is adequately described by GR. 
The parameterization of eq. (j27p applies only to inter- 
mediate length scales between the horizon scale and the 
smaller transition scale a^+^/^Y^oT. However, it implies 
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FIG. 7: The ratio of the density transfer function to the 
predicted GR case, plotted at a = 1 versus dimensionless 
wavenumber ky/al for two values of a\/a2 and two values 
of s. For a — 0.5 the transfer functions are almost identical, 
except that \D/ D{GR) — 1| is reduced by approximately 13% 
for qi/q2 — 0.5 and 30% for 01/0:2 = 2. 



that for long wavelengths and at early times, gravity re- 
duces to GR (with constant gravitational couphng). This 
assumption can be relaxed at the cost of introducing ad- 
ditional parameters, which seems premature given the 
difficulty of measuring any post-Newtonian parameters. 



Also, for wavelengths short compared with a 
and a^+^'/^v^ but large compared with the nonlinear 
scale, the gravitational couplings are constant but differ 
from GR, e.g., 7 = i for f{R). 

From the perspective of model testing, scale-dependent 
modified gravity is much more complicated than the 
scale- independent case considered in Sect. [Ill The mod- 
els have four dimensional parameters plus an exponent 
giving the time dependence. However, the situation is 
not so bleak, because structure formation depends only 
on G$(A;, t) and not on ^(k, t). In particular, matter den- 
sity perturbations on scales larger than the Jeans length 
and smaller than the Hubble length follow from integra- 
tion of 



5 + HS =^ 47rG<i,(fc, t)a^pmS 



(28) 



At early times, G* — > G and S evolves as in the GR solu- 
tion until the scale-dependent terms in eq. (|27p become 
important. The density transfer function D{k, t) given 
by eq. ([5]) is now scale-dependent at late times, implying 
a change in the shape of the matter power spectrum. It 
is easy to see that the transfer function can depend only 
on the dimensionless variables (fc^^aT, ai/a2, a). 

The most interesting new feature of scale-dependent 
modified gravity is the change of shape of the matter 
density transfer function. Figure [7] shows D(k, o = 1) 
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normalized to the GR result, obtained by numerically 
integrating eq. (|28p with initial conditions given by the 
GR result Ti?D{a) = | for a matter-dominated 

universe at a = 0.03. As expected, at large length 
scales {ky/oi ^ 1) the results converge to the GR limit. 
At short length scales, gravity is weaker than GR if 
ai/a2 < 1, leading to reduced growth; the growth is 
enhanced for ai/a2 > 1. Thus, scale-dependent gravity 
changes the shape of the matter power spectrum [43 • Ul- 
timately, measuring this scale-dependence (and doing so 
at several redshifts) can constrain scale-dependent mod- 
ified gravity theories. However, the interpretation of the 
galaxy power spectrum shape is complicated by scale- 
dependent biased galaxy formation and by the dark- 
matter-dependent transfer function (e.g., the neutrino 
fraction). Thus, while the linear growth of structure of- 
fers a potentially powerful test of GR versus modified 
gravity, it must be combined with other tests. 

The second function characterizing scale-dependent 
modified gravity, 7(fc,t), is (in our analysis, which as- 
sumes no particular Lagrangian) unrelated to G'$(fc,i). 
This function is best constrained by combining weak 
gravitational lensing and galaxy clustering measurements 
made at the same redshift. Care is required because the 
lensing amplitude is proportional to (1 + 7)<I> while the 
galaxy density is proportional to D and also depends on 
biasing. Galaxy peculiar velocity measurements could be 
used, in principle, to reduce or ideally eliminate the de- 
pendence on biasing However, one must be careful 
not to assume the velocity-density relation obtained in 
GR. The continuity equation gives 



ikdlnD 
9 In a 



HS 



(29) 



where the logarithmic growth rate dlnD/dlna is now 
scale-dependent, as shown in Figure m As in the case of 
galaxy clustering, measurement of this effect is contin- 
gent upon knowing the composition of dark matter (hot 
dark matter has a free-streaming scale, and its abundance 
determines the suppression of growth at small scales) and 
correcting for any velocity bias. 

The greater freedom allowed by scale-dependent mod- 
ified gravity models, and the fact that astrophysics (bi- 
ased galaxy formation and dark matter dynamics) may 
also introduce scale-dependence into transfer functions, 
makes it challenging to test GR using growth of struc- 
ture and weak gravitational lensing. It is likely that a 
combination of galaxy clustering, peculiar velocities, and 
weak lensing will be needed to obtain strong constraints 
on scale-dependent modified gravity theories. 



VI. MODIFIED GRAVITY VERSUS SHEAR 
STRESS 

A difference between the two longitudinal potentials $ 
and need not signal modified gravity; it might arise 
from shear stress [l3,ll3- For scalar mode fluctuations. 
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FIG. 8: Scale-dependence of the logarithmic density pertur- 
bation growth rate din D/ din a at a = 1 for scale-dependent 
modified gravity models. 



the shear stress is fully characterized by a scalar potential 
TT, such that the spatial stress tensor components are 



rpl 



p6' 



p) V'V, 



(30) 



where {p-\-p) is the background enthalpy and A = V^V^. 
In linearized GR, one of the Einstein field equations 
yields 



* - $ = l2nGa^(p + p)tt 



(31) 



All of the results obtained in Sections HTl and Hill for mod- 
ified gravity apply equally to GR with shear stress if 7 is 
replaced by 127rG'a^(p-|-p)7r/$. 

In standard cosmology, the only significant source of 
shear stress is relativistic neutrinos after neutrino de- 
coupling during the radiation-dominated era. For long 
wavelengths during the radiation-dominated era, neu- 
trino shear stress gives [13] 



5 V P + P 



(32) 



During the matter-dominated era, 7 — 1 cx a^^ and shear 
stress is unimportant at late times in the ACDM model. 
It is also unimportant in simple quintessence models be- 
cause shear stress vanishes for linear perturbations of a 
minimally-coupled scalar field. 

Shear stress might nonetheless be important if cosmic 
acceleration is driven by an imperfect fluid. Without 
specifying the dynamics of this fluid, few constraints can 
be placed on tt. One possible bound comes from the 
dominant energy condition, which states that each of the 
eigenvalues of T*^ must be smaller in absolute value than 
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the energy density. If this condition holds, then eqs. (1301) 
and (|3ip can be combined to give rather weak bounds on 
^-Z^- 1. 

Additional constraints follow from the initial-value 
constraints of GR and energy-momentum conservation, 
which for a spatially flat background become [l^ 



5 + ma 

u + (1 - ■icl)nu 



AwGa'^ {p + p) {5 + mu) , (33a) 

AT:Ga^{p + p)u , (33b) 

3* - k^u , (33c) 

c^J + o- + $ - fcV . (33d) 



The first of these equations is the usual Poisson equation 
in conformal Newtonian gauge. The density and veloc- 
ity potential perturbations are defined from the energy- 
momentum tensor components by T^q = —p — {p + p)S, 
r'^j = —(p + p)'S/iU, while the entropy perturbation is 
defined by cr = {Sp — cl5p)/{p -I- p) with = dp/dp. 
For a single perfect fluid like cold dark matter before its 
trajectories intersect, u = 0. However, in general cr 7^ 
for a multi-component fluid, e.g. dark matter and a non- 
constant dark energy. 

The freedom introduced by entropy and shear stress 
perturbations is, unfortunately, sufficient in principle to 
reproduce any observations of large-scale structure and 
gravitational lensing. Consider, for example, perfect 
measurements of galaxy peculiar velocities everywhere 
and at all times assuming that galaxies exactly trace cold 
dark matter. Then, eq. (|33d[) with = ct = vr = for 
CDM yields $(k, i). Assume furthermore that complete 
and ideal gravitational lensing measurements are avail- 
able to yield $(k, t) + ^(k, t). Now, the GR initial- value 
constraints p3ap - (|33bp suffice to yield ^(k, t) and u(k, t) 
for the multi-component fluid of dark matter and dark 
energy. Requiring this fluid to obey energy-momentum 
conservation (j33cp - (|33dp yields a and tt. In short, per- 
fect measurements of $ and can be used to deter- 
mine a and tt of the combined fluid of dark matter and 
dark energy. When combined with measurements of the 
dark matter density and velocity fields (assuming galax- 
ies trace dark matter), one can, in principle, determine 
the energy-momentum tensor components for dark en- 
ergy. One can think of this energy-momentum tensor for 
dark energy as the difference between the Einstein ten- 
sor for the observed metric and the energy-momentum 
tensor of the ordinary matter 11 1 . 

This approach can be used, in principle, to determine 
the dark energy entropy and shear stress needed to ex- 
plain any observations of large-scale structure (including 
peculiar velocities) and weak lensing — even if there is no 
dark energy, but instead gravity is modified. In effect, the 
two observable metric fields $ and 5* can be exchanged 
for either cr and tt (GR with exotic dark energy) or G$ 
and 7 (modified gravity). 

Although the initial-value constraints of GR can be 
used to determine the properties of dark energy, they can- 
not be used to model modified gravity. As we have seen 
in previous sections, modified gravity leads generically to 



a modified Poisson equation with a variable gravitational 
coupling G$(/e, t). In order not to break local Lorentz in- 
variance by the selection of a preferred frame, the other 
components of the Einstein equation should also be mod- 
ified. For example, in f{R) theories the left-hand side of 
(|33bp is multiplied by (1 4- Jb) while the right-hand side 
acquires an extra term 5/7?,. Neglecting these modifica- 
tions leads to violations of eq. ^ on large scales. For 
example, Caldwell et al. [l3] modeled the dynamics of 
modified gravity with three recipes that assume the va- 
lidity of some combination of eqs. (j33a[) and (j33bp . While 
recipes Rl and R2 satisfy the conservation of C,, recipe R3 
does not. As a result, it leads to a different prediction for 
the gravitational potentials and therefore the ISW effect. 



VII. DISCUSSION 

Parameterizing modified gravity theories in cosmol- 
ogy is much more difficult than parameterizing post- 
Newtonian gravity in the solar system because the gravi- 
tational potentials and ^ in the GR limit are not static 
Coulomb potentials in cosmology. In GR, on scales larger 
than the Jeans or nonlinear length scales, = <i> factor- 
izes into a product of a time-dependent growth function 
and a spatially-varying curvature perturbation. If this 
factorization persists in modified gravity theories, then 
a simple post-Newtonian parameterization can be ob- 
tained. This is the approach we introduced in Section 
im It has the practical virtue of yielding easily calcu- 
lated predictions for all observables in the linear regime, 
including the growth of matter clustering, peculiar veloc- 
ities, microwave background anisotropy, and weak gravi- 
tational lensing. 

Introducing the Eddington parameterization ^E" = 7$, 
where 7 depends on time but not on space, we showed 
that structure grows faster (gravity is stronger) in models 
with 7 < 1. However, the shape of the transfer functions 
(i.e., their dependence on spatial wavenumber) for these 
scale-independent models is unchanged compared with 
GR. 

Unfortunately, realistic modified gravity theories such 
as the f{R) models have scale-dependent effects and can 
no longer be described by only one post-Newtonian quan- 
tity 7. Instead, the strength of gravity, described by 
G$ in the Poisson equation, can vary independently of 
7. Even so, because galaxy clustering and dynamics de- 
pends on G$ but not on 7, while weak lensing depends 
on 7, observations could, in principle, measure modified 
gravity parameters assuming there is no dark energy. 

If dark energy is complex enough to require two addi- 
tional fields to characterize its stress tensor (e.g., shear 
stress potential and entropy), then it appears that there 
is enough information in the dark energy model to ac- 
count for any $ and ^' without modifying gravity. One 
cannot prove gravity is modified unless one can account 
for all significant contributions to the stress-energy ten- 
sor. 
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Thus, our hope to describe all modified gravity mod- 
els with two parameters, yielding predictions measurably 
different from all dark energy models, has not been real- 
ized. Distinguishing modified gravity from dark energy 
will require making additional assumptions. 

Nevertheless, the s) parameterization of scale- 
independent modified gravity presented in Section|TTl and 
the {ai, a2, (3i^ P2) parameterization of scale-dependent 
modified gravity models presented in Section |Vl are still 
useful for characterizing observational data. If measure- 
ments of galaxy clustering, peculiar velocities, and weak 
lensing are all consistent with f3 — and ai — a2 — Pi — 
(32 = 0, for example, then modified gravity and exotic 
dark energy models can both be excluded. If measure- 
ments require nonzero parameters, however, dark energy 
and modified gravity remain viable explanations until ad- 
ditional assumptions are made to distinguish them, e.g. 
restriction of the Lagrangian to a particular form. 

A generic prediction of modified gravity theories in cos- 
mology is that the gravitational coupling in the Pois- 
son equation should vary with time and with length scale. 



Departures from GR could be important not only in the 
linear regime of cosmological perturbations but perhaps 
also in the nonlinear regime (albeit on scales much larger 
than the solar system). Nonlinear effects may allow mod- 
ified gravity to be distinguished from exotic dark energy, 
assuming that the dark energy fluctuations are small. For 
this reason it would be valuable to perform N-body sim- 
ulations of structure formation using variable G$, ex- 
tending previous work [4^ to the scale-independent and 
scale-dependent modified gravity models discussed in the 
current paper. 
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